Cooling and squeezing the fluctuations of a nanomechanical beam by indirect 
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We study cooling and squeezing the fluctuations ol a nanomechanical beam using quantum feed- 
back control. In our model, the nanomechanical beam is coupled to a transmission hne resonator via 
a superconducting quantum interference device (SQUID). The leakage of the electromagnetic field 
from the transmission line resonator is measured using homodyne detection. This measured signal 
is then used to design a quantum-feedback-control signal to drive the electromagnetic field in the 
transmission line resonator. Although the control is imposed on the transmission line resonator, this 
quantum-feedback-control signal indirectly affects the thermal motion of the nanomechanical beam 
via the inductive beam-resonator coupling, making it possible to cool and squeeze the fluctuations 
of the beam, allowing it to approach the standard quantum limit. 

PACS numbers: 85.25-j, 03.65.Ta, 42.50.Lc 
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I. INTRODUCTION 

Nanomechanical oscillators have recently attracted 
considerable attention for their possible applications in 
quantum informat_ion_and c[uantum measurement (see. 
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cillator is also a promising device for studying macro- 
scopic quantum effects in mechanical systems (see, e.g., 
Refs. [1 0, i i, S, H, S i, i, [S)- Using current ex- 
perimental techniques (see, e.g., Refs. [1, Ull, El]): high- 
frequency nanomechanical oscillators (a;/27r ~ 1 GHz) 
with quality factors Q in the range of lO'^-lO^ can be re- 
alized at low temperatures T on the order of mK. When 
the vibrational energy fku of the nanomechanical oscilla- 
tor becomes smaller than the thermal energy ksT, the 
oscillator can be said to work in the quantum regime. 

To observe quantum behavior in nanomechanical os- 
cillators, e.g., quantum fluctuations or squeezing effects, 
the oscillator must be cooled to extremely low tempera- 
tures to approach the standard quantum limit. There 
have been numerous studies, both theoretical and ex- 
perimental (see , e.g., Refs. lie 

M MM Sip, '25', 'Seljp? 

M m [sa [3H33I39, 40, 4l|)^ investigating the cooHng 
of the fluctuations of nanomechanical oscillators. Many 
of these studies focus on optomechanical systems (see, 

e.g., Refs. [H, 0, m M M, MMMMMm), 

where an oscillating cantilever or an oscillating micro- 
mirror is modelled as a harmonic oscillator. There are 
two approaches in optomechanical cooling: passive cool- 
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la [M [13, H [H, Hi and active cool- 
23l |. In passive cooling techniques, the 
mechanical oscillator is self-cooled by the dynamical 
back-action, e-g., the radiation-pressure-induced back- 
action di, [H, [m m [n m 13 

coming from the mirror 
surface of the optical cavity. In fact, for a high-finesse 
cavity, the photons reflected from the mirror of the cav- 
ity transfer momentum and induce additional damping 
to the mechanical oscillator. In active cooling techniques, 
the reflected signal coming from the mechanical oscillator 
is sent to an electronic circuit, e.g., a derivative circuit, to 
provide a modulating signal, which is then used to control 
the back-action force imposed on the mechanical oscilla- 
tor. Since the cooling effect can be actively controlled by 
tuning the feedback gain obtained in the control circuit, 
this is called an active cooling strategy. 

Although it has recently been reported that ground- 
state cooling [H, [2^, [2^, [23] could be realized in optome- 
chanical systems, it is difficult to observe the macroscopic 
quantum effects of the mechanical oscillators in these op- 
tomechanical systems using current experimental condi- 
tions. The main difficulty comes from the fact that the 
characteristic oscillating frequency of the mechanical os- 
cillator in these systems is not high enough (typically on 
the order of kHz or MHz), and the corresponding effective 
temperature to observe the quantum effects is extremely 
low (typically on the order of nK or /iK), which is difficult 
to realize in present-day experiments. 

Besides optomechanical cooling, a nanomechanical os- 
cillator can also be embedded in an electronic circuit and 
cooled by coupling it to an electronic system [2^, [2^, 
Possible strategies include nanomechanical oscilla- 
tors coupled to superconducting single-electron transis- 
tors [31, 32], quantum dots 1,33. .34L Josephson-junction 
superconducting circuits 135], l36l 1371 . l38l l39l | , or transmis- 
sion line resonators (40l . l4l| . Compared with mechanical 
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oscillators in optical systems, a high-frequency oscillator 
can be realized more easily in electronic systems. Indeed, 
it has been reported that nanomechanical beams [HI, 
with frequencies in the regime of GHz have been real- 
ized, and these beams seem to be suitable for integration 
in an electronic circuit. Since the effective temperature 
of such a mechanical oscillator can be in the mK regime, 
it should be possible to observe quantum behavior in this 
case. 

Like optomechanical systems, active feedback controls 
can be introduced to cool the motions of the nanome- 
chanical oscillators in electronic systems. In the theo- 
retical proposal in Ref. [33 |. a nanomechanical resonator 
is capacitively coupled to a single-electron transistor to 
measure the position of the resonator. The information 
obtained by the quantum measurement is fed into a feed- 
back circuit to obtain an output control signal, which 
is then imposed on a feedback electrode to control the 
motion of the resonator. In this strategy, the quantum 
measurement and the designed feedback control intro- 
duce additional damping effects on the resonator, which 
are helpful for cooling the motion of the resonator. 

In a recent experiment [4^ . a nanomechanical beam 
acted as one side of a superconducting quantum interfer- 
ence device (SQUID), and the voltage across the SQUID 
was measured which can be used to detect the motion of 
the beam. Motivated by this experiment, here we study 
the coupling between such a system [4^ and a trans- 
mission line resonator. A single-mode quantized elec- 
tromagnetic field provided by this resonator [431 could 
be detected by a homodyne measurement [44| . The 
quantization of this coupled beam-SQUID-resonator sys- 
tem has been addressed in the literature (see, e.g., 
Rcf. [4^, [4^, 1131), and theoretical analysis shows that 
such a device can be used to detect the motion of the 
beam [l^]. Here, we would like to concentrate on a dif- 
ferent problem: how to design a quantum feedback con- 
trol from the output signal of the homodyne detection to 
drive the motion of the beam? Different from previous 
work, such as the one in Ref. [3^ . the quantum feedback 
control proposed here is imposed on the transmission line 
resonator, not on the beam, and indirectly controls the 
motion of the nanomechanical beam via the coupling be- 
tween the transmission line resonator and the beam. By 
adiabatically eliminating the degrees of freedom of the 
SQUID and the transmission line resonator, the designed 
feedback control could introduce anharmonic terms in the 
effective Hamiltonian and additional damping terms for 
the nanomechanical beam, leading to both cooling and 
squeezing of the fluctuations of the beam. 

This paper is organized as follows: in Sec. [TTl we 
present a coupled system composed of: (i) a transmis- 
sion line resonator, (ii) an rf-SQUID, and (iii) a nanome- 
chanical beam. The open quantum system model and the 
quantum weak measurement approach used to design the 
feedback control are investigated in Sec. IIIII The quan- 
tum feedback control design is presented in Sec. II VI and 
our main results about squeezing and cooling the fluctua- 



tions of the beam are discussed in Sec.|Vl The conclusion 
and discussion of possible future work are presented in 
SedVH 



II. SYSTEM MODEL AND HAMILTONIAN 

We mainly focus on a physical model in which a 
doubly-clamped nanomechanical beam, a rf-SQUID, and 
a transmission line resonator are inductively coupled (see, 
e.g., Ref. [45., ) . The quantum electromechanical 

circuit and the corresponding equivalent schematic dia- 
gram are shown in Fig. [T] In this circuit, the mechanical 
oscillator, i.e., the clamped nanomechanical beam, is in- 
tegrated into the rf-SQUID with a Josephson junction 
having a critical current Ic and a capacitance C . Here, 
the displacement of the beam in the plane of the loop 
with a small amplitude x around its equilibrium position 
changes the area of the loop, and thus influences the total 
magnetic flux $ threading the loop. There is an applied 
external flux threading the loop. The rf-SQUID inter- 
acts with a nearby transmission line resonator (TLR) , via 
their mutual inductance. The additional magnetic flux 
provided by the quantized current in the transmission 
line resonator is 



add 



= $T(-«a + ia^)i 



where $t is a constant; a and are the annihilation and 
creation operators of the quantized electromagnetic field 
in the transmission line resonator. Here, we ignore the 
small change of <I>t caused by the oscillation of the beam. 
For this superconducting circuit, the total magnetic flux 
<I> threading the loop of the rf-SQUID is given by: 



$ = $e + Blx + ^T{-ia 



LI, 



where L and / are, respectively, the self-inductance and 
the current in the loop; / is the effective length of the 
beam; and B is the magnetic field threading the loop 
of the rf-SQUID at the location of the nanomechanical 
beam, and is assumed to be constant in the region where 
the beam oscillates. 

The total Hamiltonian of this coupled electromechan- 
ical system can be written as [i^: 



H = Ho + Hi 



(1) 



with the Hamiltonians: 



Ho 



P 

2m 



hujTO^a + hu{t)(a^ + a), (2) 



2ttBI 

— K[—ia + la' ) — X 



+2Uo(3LCOS(t)+ 



(3) 



where $0 is the flux quantum; lvm is the oscillating fre- 
quency of the doubly-clamped beam; and 4>e are related 
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FIG. 1: (color online) (a) Schematic diagram of a transmission 
line resonator (TLR), in blue, and a SQUID-nanomechanical 
beam system, in beige and dark blue. Here, "JJ" represents 
a Josephson junction, and Ei-n, Eout are, respectively, the in- 
put and output electromagnetic fields of the transmission line 
resonator, (b) Diagram of the quantum circuit. 



to the normalized total flux and external flux: 



= 27r 

\$o 2 



' $0 2 



The normalized system parameters C/q, /3l and n in 
Eq. Ill) are given by: 



3_ 

87ri' 



2ttLIc 



The observables p and q in Eqs. ([2]) and ([3]) are, respec- 
tively, the conjugate observables of x and <& represent- 
ing the momentum of the beam and the charge on the 
Josephson junction. The term hu{t){a^ + a) in Eq. ^ 
is an interaction Hamiltonian between the transmission 
line resonator and the external control field, where the 
time-dependent function u(t) can be designed according 
to the desired goal. 



When /3l > 1 and 
2TrBlx 



-|- K{—ia + ia^) 



«1, 



(4) 



the Hamiltonian Hi represents a double-well potential 
near = 0, and the two lowest eigenstates, \L) and \R), 
correspond to two current states with opposite circulat- 
ing currents in the loop of the rf-SQUID, which are far 
separated from higher-energy eigenstates. At sufficiently 
low temperatures, only the two lowest eigenstates \L), 
\R) contribute. Thus, the rf-SQUID can be modelled as 
a two-level system, and the Hamiltonian Hi of the rf- 
SQUID can be re-expressed as [46| : 

he / 2TrBl , . . .. \ ^ HA^ 

Hl = — {<Pe + X + K{-ia + ia') j C7z 

where ax and az are the x-axis and z-axis Pauli operators 
in the basis of \L) and \R); e, A are real parameters that 
determine the energy difference between the two minima 
of the double- well potential and the tunnelling amplitude 
between the wells, respectively. Under the condition that 

< - 1) « 1, 

then e and A can be approximately given by [isj : 



V6(/3l - 1), 



A = 3[/o ( 1 - 

PL 



(5) 



Letting the external flux 0e = 0, we can rewrite the 
Hamiltonian Hi in the qubit basis as: 



Hi 



hujs 



hneBl 



+ ——{-ia + ia))ax, (6) 



where ws = A, and Ux, <Jz are the corresponding x-axis 
and z-axis Pauli operators in the qubit basis: 



\R). 



Here, we assume that the oscillation frequency of the 
beam is high enough such that lom, ojs, and ujt are com- 
parable. Then, under the rotating-wave approximation, 
and with Eq. ([6|), the total Hamiltonian in Eq. ([T]) be- 
comes: 



H = 



hujs 



+ hujMb^b + hujTO^ a 



+hu{t){a^ + a) + hgMs{b<y+ + <J-b'') 



(7) 



where the coupling strength gus between the mechanical 
oscillator and the rf-SQUID is: 



9ms — ffMcch-SQUID 



-KtBl 



4 



and the coupling strength gsx between the rf-SQUID and 
the transmission hnc resonator is given by: 



where X^^h' ^baih and X^^^^^^ are, respectively, the en- 
vironmental operators interacting with the transmission- 
line resonator, the beam and the rf-SQUID; and Mr 
are two constants which determine the dcphasing and re- 
laxation rates of the rf-SQUID. Furthermore, let us con- 
sider a bosonic model of the environment, and assume 
that the interactions between the system degrees of free- 
dom and the environmental degrees of freedom are linear 
interactions. Then, under the rotating- wave approxima- 
tion and the Markovian approximation, we can obtain the 
following quantum stochastic differential equation for a 
system observable X (see Appendix[X]for the derivation): 



.(2) 



(3) 



9ST — ffSQUID-TLR 



2h' 



The annihilation and creation operators b and of 
the fundamental oscillating mode of the nanomechanical 
beam are defined by: 



mujM 

2h ' 



1 



\/2hmujM 
1 

^ — 

\/2hmujM 



P, 



p. 



Furthermore, let us assume that the frequencies of the 
rf-SQUID, the beam, and the transmission line resonator 
satisfy the conditions: 

gMS ^ Aj\/5 — LOs^l^M, gsT <C AsT = UJS — ^T- (8) 

Then, in this large-detuning regime [HI, the follow- 
ing transformation can be introduced to diagonalize the 
Hamiltonian H in Eq. ([7]): 



U 



exp 



gMS 



A 



MS 



^ST 



In fact, under the condition given in Eq. ([8]), we can 
obtain an effective Hamiltonian: 

~ hujMb^b + hujTCL^ a + hu{t){a'^ + a) 

+ ^a. + h(4^b^b+^a^a) a. 
2 \Ams Ast J 

f- ( gMSgST , gMS9ST\ , t , -.t n 

by expanding UHW to first order in gMS / Ams a-nd 
gsr/AsT- 



III. INTERACTION BETWEEN THE SYSTEM 
AND ITS ENVIRONMENT 

A real physical system inevitably interacts with the 
external degrees of freedom in the environment. Such in- 
teractions introduce noise to the system. There are three 
kinds of noise that should be considered here: the ther- 
mal noises on the nanomechanical beam and the trans- 
mission line resonator, as well as the electromagnetic fluc- 
tuations on the rf-SQUID caused by the nearby electro- 
magnetic elements. 

The interaction Hamiltonians between the transmis- 
sion line resonator, the beam, the rf-SQUID and their 
environments can be described as: 



rr(l) _ 



H 



(2) _ 



bath' 
(2) 



H.. 



int 

(3) 



(a 4 



dX = 



+ ^(nAf + l)(6t[A,fe] + [6t,A]&) 
+ l^nM{b[X,b^] + %Xp) 
+lL{a^[X,a] + [a\X]a) 



/^dAl[X, 



/Tj^[a\ X]dAin, (9) 



where 7m, 7s, 7t are the damping rates of the mechan- 
ical beam, the rf-SQUID and the transmission line res- 
onator under the Markovian approximation; 



1 



(10) 



is the average photon number of the beam in thermal 
equilibrium with the environment at temperature T. To 
simplify our discussions, when Eq. Q was derived, we ne- 
glected environment-induced thermal excitations on the 
transmission line resonator and the rf-SQUID (these ex- 
citations could indeed be neglected with the parameters 
given in Eqs. (EH) and ^ in Sec. |V|. The leakage of 
the transmission line resonator could be detected using 
a homodyne detection with detection efficiency 77, where 
dAin represents a quantum Wiener noise 50] satisfying: 



dAidA;,, = dA„dAi,, = dAidA;., = 



dA;„dA 



dt. 



Here, we only keep the fluctuation terms caused by 
the measurement and average over the other fluctua- 
tions, because the evolution of the coupled beam-SQUID- 
resonator system is conditioned on the measurement out- 
put, which depends on the measurement-induced fluctu- 
ations. The corresponding measurement output of the 
homodyne detection can be expressed as i47i] : 



dYt = 



(11) 



Note that this measurement output depends on the input 
noise and the electromagnetic field of the transmission 
line resonator. 
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IV. QUANTUM FILTERING AND QUANTUM 
FEEDBACK CONTROL 



There are two possible ways to design a quantum feed- 
back control protocol |5ll, H^] based on the measurement 
output. One approach is to directly feed back the out- 
put signal to design the quantum feedback control signal, 
which leads to the Markovian quantum feedback con- 
trol [5^. Another approach, which is called quantum 
Bayesian feedback control [54|, can be divided into two 
steps: the first step is to find a so-called quantum filter- 
ing equation [sl, IH, [s^l to give an estimate of the state 
of the system from the measurement output; the sec- 
ond step is to design a feedback control signal based on 
the estimated state. The possibility for a "control prob- 
lem" to be divided into these two steps, i.e., a separate 
filtering step and a control step, is called the separation 
principle in control theory, which has recently been devel- 
oped for quantum control systems [56]. Compared with 
the Markovian quantum feedback control, the quantum 
Bayesian feedback control can be applied to more gen- 
eral systems. In our proposal, we will design the control 
using Bayesian feedback control. 



Based on quantum filtering theory, which has been well 
developed in the literatures [55, 56|, we can obtain the 
following stochastic master equation for the estimated 
state p: 



dp 



-j\HcS, p]dt + -fMnuT) \b'^] pdt 

+1M {fiM + 1) V[b]pdt + jsMlV[a,]pdt 
+-fsM^V[a_]pdt + jT'D[a]pdt 
+Vvi^'H[a]p {dYt - VTf¥{{a^ + a))dt) , (12) 



where 



p = E{p\yt) 



(13) 



is defined as the conditional expectation of the density 
operator p for the coupled beam-SQUID-resonator sys- 
tem under the von Neumann algebra [56| : 



yt=vN{Ys\h<s<t}, 



(14) 



spanned by the measurement outputs; (A) — tr(Ap) is 
the average of A under p; and the superoperators 'D[c]p 
and 7i[c]/5 are defined by: 



V[c]p 

n[c]p 

The increment 

dW = dYt 



cpc' 



- t - 
-c' cp - 



~ t 
-pc'c. 



cp + pc' 



{{c + c^))p. 



/rijT{{a^ + a))dt 



(15) 



in Eq. ()12p is the innovation updated by the quantum 
measurement, which has been proved to be a classical 
Wiener increment satisfying: 



E{dW) = 0, {dWf 



dt. 



for homodyne detection (see, e.g., Ref. |56|). 

The von Neumann algebra J^t defined in Eq. rep- 
resents the information obtained by the quantum mea- 
surement up to time t; thus the conditional expectation 
p defined in Eq. (fTS]) is the best estimate of the system's 
state obtained from the measurement output. 

Generally, the stochastic master equation, i.e., the fil- 
ter equation, is difficult to solve. However, for the sys- 
tem we discuss here, the stochastic master equation (fT2)) 
is equivalent to a set of closed equations under the semi- 
classical approximation (see Eqs. (jB3p and (jB4p in Ap- 
pendix [B|). This set of equations can be integrated by a 
Data Acquisition Processor (DAP) (see, e.g., Ref. [51]), 
which is composed of a Digital Signal Processor (DSP) 
and analog/digital and digital/analog signal converters. 
Such a Data Acquisition Processor works as an integral 
estimator of the dynamics of the system state and gives 
the output signals (xt) and (pr)- The output signals of 
the estimator are fed into a feedback controller (a linear 
amplification element) to obtain the following feedback 
control signal: 



u{t) = V {-Vx{xt) + Vpipr)) 



(16) 



where and Vp are the feedback control gains that can 
be chosen according to the desired goal, and xt, Pt are 
the normalized position and momentum operators of the 
transmission line resonator defined by: 



XT 



PT 



(17) 



We replace u(t) in Eq. ^ and Eq. ^ by Eq. ^ to 
obtain new dynamical equations. In this case, we indeed 
control, simultaneously, the evolutions of the transmis- 
sion line resonator and the estimator. The schematic di- 
agram of the feedback control circuit is shown in Fig. [21 
From the definition of p, it can be shown that the 
control of the coupled system given by Eq. Q is equiv- 
alent to the control of the estimator given by Eq. (|12p . 
Thus, in the following discussion, we will focus on how 
to control the quantum filtering equation (|12|) . 

If the damping rates 75 and 7t of the rf-SQUID and 
the transmission line resonator are large enough such that 



7S, 7T > IMUM, 



(18) 



we can adiabatically eliminate [58|, |59| the degrees of free- 
dom of the rf-SQUID and the transmission line resonator 
to obtain the following reduced stochastic master equa- 
tion and the measurement output for the nanomechanical 
beam (see Appendix [Bl for the derivation): 



dpM = 



dYt 



PM]dt + "fMnAiT^ib^PMdt 

+lM{nM + l)V[b]pMdt 
+jTT>[Cib + C2b^]pMdt 
+Vm^'H[Cib + C2b'<]pMdW, 
^/^{{axb + alb^))dt + dW, (19) 
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_L_ 

rf-SQUID 



Nanomechanical beam 



FIG. 2: (color online) Schematic diagram of the feedback 
control circuit. The "DSP" denotes a Digital Signal Proces- 
sor (DSP) which works as the integral estimator of the sys- 
tem state by solving Eqs. (fB3)) and (|B4)) . The "A/D" and 
"D/A" represent the analog/digital and digital/analog signal 
converters. The output of the estimator is fed into a linear 
amplifier circuit to obtain a control signal which is further 
used to drive the input electromagnetic field of the transmis- 
sion line resonator and the input of the estimator. 



Equation (jl9p shows that the quantum measurement 
and feedback control introduce extra damping and fluc- 
tuation terms for the beam (the third and fourth lines 
in Eq. (|19p ). These damping terms are important for 
squeezing and cooling the fluctuations of the beam. 

V. SQUEEZING AND COOLING THE 
FLUCTUATIONS OF THE NANOMECHANICAL 
BEAM 

In order to study the squeezing and cooling effects on 
the nanomechanical beam induced by the quantum feed- 
back control, let us first define the normalized position 
and momentum operators of the nanomechanical beam: 

XM ^ ^(b + b^), PM ^ ^i~ib + ib^). 

Then, from the reduced stochastic master equation (|19p . 
we can study the evolutions and the corresponding sta- 
tionary values of the variances 

= (4f) - (^M)^ V^pM = {pIi) - {PM? (22) 
of xm and pm- 

A. Squeezing 



where the reduced Hamiltonian H^^^'^ is given by: 

= huMb'^b + h^Mb' + hCMb''^ 

+hu{t){a^h + alh'^), (20) 

and the reduced effective control on the beam is: 

u[t) = -Vx{[axb + a^.b'^)) + Vp{{aph + ayj^)) . 

The parameters a^, ap, tM can be expressed as: 

Olx = G1+C2, 

ap = —iCi + 16*2 , 
Cm = OJTC2C1 — igMTC2, 



where gMT is given by Eq. (jB2p and Ci, C2 are given by: 





) - 


9mt 
X 


[(» 




) = 


9MT ^ 

X 


Vp - 




) = 


IT , 

— {IT + 



As shown in Eq. ([20]) . there is a two-photon term 

^S.Mb'^ +fi£,*Mb'^ ^ in the effective Hamiltonian h'^^^ , which 
leads to squeezing in the fluctuations of the beam. With- 
out the quantum feedback control, i.e., Vx = Vp = 0, S^m 
would be zero and the two-photon term vanishes. 



The nanomechanical beam can be described by the 
conjugate variables xm and pu, and is in a squeezed 
state if the corresponding variances of these variables de- 
fined in Eq. satisfy 14^ < ^■72 or Vp^^ < h/2, i.e., 
the variance of one of the two conjugate variables is below 
the standard quantum limit (see, e.g., Refs. [iSllllIll)- 
Owing to the uncertainty principle which requires that 
^Af^M — ^^/4, squeezing the fluctuations of one of the 
two conjugate variables would lead to the dispersion of 
the fluctuations of the other conjugate variable. 

If we choose the feedback control gains Vx and Vp to 
satisfy the conditions 



It 

< < 1, 

7T + 4:Vp 

LOT — 2Wt 

< — < 1, 

LOT 

lom 



(23) 



then the stationary variances V^^^ and Vp^^ can be esti- 
mated as: 




PM 



h 1 

2v!^' 



(24) 



where 



lom - 2ReC_ 



M 



lom 



2ReCM' 



Re^Af 



^TgMT 

x' 



{vl-vl). (25) 
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It is shown in Eq. ([24)1 that the parameter ^ determines In this case, we have Vp > v^. Then, from Eqs. p4|) 

and (1^ . it can be verified that VT < Vlf , which co- 
incides with the numerical results in Fig. Wia) (the blue 
dashed line for V^^ Jh is below the green triangular line 



the tradeoff of the squeezing effects between V^m and 
■ When ^ > 1, the fluctuation of the momentum pM 
of the nanomechanical beam is squeezed. However, when 
^ < 1, the fluctuation of the position xm of the nanome- 
chanical beam is squeezed. The parameter 77, i.e., the 
measurement efficiency of the homodyne detection, de- 
termines the minimum uncertainty that can be reached. 
For a quantum weak measurement with a high efficiency 
77 such that 



< 77 < 1, 



(2nM + 1)^ 



where um is the thermal excitation number of the beam 
given in Eq. pi7|) . it can be verified that 

— < « — 

4 '^XM^PM 4^ 



where V^^^ and T^"^ are the "uncontrolled" stationary 
variances that are obtained from Eq. (|19|) by letting = 
Vp = 0. It should be pointed out that the product of the 
uncertainties of xm and pm given by Eq. \24^, i.e., 



"xm^'pm 4^^' 

corresponds to the Heisenberg uncertainty limit of a quan- 
tum system under imperfect quantum weak measurements 
(see, e.g., Ref. [63j). When the measurement efficiency 
77 tends to unity, the traditional Heisenberg uncertainty 
limit /i^/4 is recovered. 

To show the validity of our strategy, let us show some 
numerical examples. The system parameters are chosen 
as gi]: 

L = 3.38 X 10-" H, C = 7.4xl0i^F, 
Ic = 10 fiA, m = 10-16 kg, ^ = 0.6, 
WA//27r = 1 GHz, B; = lTx^m, 
Q = 104, T = 100 mK, 0e = 0, 
7S jl-K = 100 MHz, 7T /27r = 20 MHz, 
gsTl2-K = 20 MHz, ut/2ti ^ 4.3 GHz. (26) 

From the above parameters, it can be calculated that 

7M/27r « 0.1 MHz, t^5/27r « 6.3 GHz 
gMs/27r « 73 MHz, gMr/'^'^ ~ 4.9 MHz. (27) 

In our numerical results, the stationary variances are cal- 
culated from the dynamic equation ((T9)) . and the feed- 
back control parameters Vj. and Vp are chosen to satisfy 
Eq. In fact, in Fig. O we choose and Vp such 

that 



for Vp^j /h). It means that the fluctuation of the position 
of the beam is squeezed. Meanwhile, the numerical re- 
sults in Fig. [3Kb) show that the variances V^^^ and Vp^^ 
under control are much smaller than the variances Vl""^ 
and Vp^j without control, which means that the designed 
feedback control reduces the variances of the position and 
momentum of the beam. The product of the variances 
under control V^,VJ:,, could even be reduced to be close 
to the Heisenberg uncertainty limit fi.^/4. 
In Fig. m we choose Vx and Vp such that 



Vx=0.5ujT, 0.3<^<0.5. 

LlJt 



(29) 



In this case, we can calculate from Eqs. (|24l) and ([251) that 
Vlf > VJr , which coincides with the numerical results 

Xm — Pm ' 

in Fig. [If a) (the green triangular line for Vp^ Jh is below 
the blue dashed line for Vlf Ih). It means that the fluctu- 

Xm I > 

ations of the momentum of the beam are squeezed. More 
interestingly, the numerical results in Fig.^^a) show that 
the variance ofpu is squeezed to be less than the standard 
quantum limit, i.e., Vp^^ < h/2. Numerical results in 
Fig.[4Kb) show that the product V^^Vp^^ of the controlled 
variances is much smaller than the product V^^Vp^ of 
the uncontrolled variances, which means that the vari- 
ances of the position and momentum of the beam could 
be reduced under the designed feedback control. Indeed, 
as shown in Fig.[3Kb), the product of the variances could 
be reduced to be close to the Heisenberg uncertainty limit 

hy4. 



B. Cooling 

Further, let us investigate the cooling of the fluctua- 
tions of the nanomechanical beam. The cooling effect 
can be estimated by the average photon number of the 
nanomechanical beam 



1 

2ft' 



(VxM +VpJ-- + ^{VixM) + ViPM)) 



1 



(30) 



with 



0.5wT, 0.5<^<1. 



UJT 



(28) 



Xm = Edwi(xM)), Pm = Edwi (pm ) ) , 
V{xm) = Edwi(xM)^) - xIj, 

V{pm) = Edw{{pM)^)-PM- 

Here, Edw means that the expectations and variances of 
(xm) and {pm} are about the Wiener noise dW. 

From the system parameters given in Eqs. ([26]) and 
([?71) , it can be verified that the controlled stationary ex- 
pectations x%j, p\,j and the classical fluctuations ^^^a/)' 
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FIG. 3: (color online) Squeezing the position-fluctuations of 
the beam with the parameters given in Eqs. (|26p . (|27p . and 
(|28l) . (a) Variances of the position and momentum of the 
beam (in units of h) versus the normalized control parame- 
ter Vp = Vp/ujT- Recall that Vp is the linear feedback gain 
(see, Eq. (|16p ). The green line with triangles and the blue 
dashed line represent, respectively, the controlled variances of 
the position V^j^,j /h and momentum V^^^ /h of the beam. The 
red line with asterisks and the blue solid line represent the 
uncontrolled variances V^^^ /h and Vp^^ /h of the beam, which 
coincide because the beam without control is in a coherent 
thermal state with equal variances of position and momen- 
tum, (b) Products of the variances (in units of It?) versus 
the normalized control parameter Vp — Vp/cur- The green 
line with asterisks and the blue solid line denote the uncon- 
trolled trajectory Vxm^pm /^^ ^^'^ controlled trajectory 
Vxm'^pm /^'^ ' respectively. The red dashed line at 1/4 repre- 
sents the Heisenberg uncertainty limit 



FIG. 4: (color online) Squeezing the momentum-fluctuations 
of the beam with the parameters given in Eqs. (|26p . (127^ . 
and (|29|) . (a) Variances of the position and momentum of 
the beam (in units of h) versus the normalized control pa- 
rameter Vp = Vp/uJT- The green line with triangles and the 
blue dashed line represent, respectively, the controlled vari- 
ances of the position V^^Jfi and momentum Vp^Jh of the 
beam. The red line with asterisks and the blue solid line rep- 
resent the uncontrolled variances V^^/h and Vp^Jfi of the 
beam, which coincide because the beam without control is in 
a coherent thermal state with equal variances of position and 
momentum, (b) Products of the variances (in units of fi^) 
versus the normalized control parameter Vp = Vp/ujT- The 
green line with asterisks and the blue solid line denote the 
uncontrolled trajectory V^'^'J /fi^ and the controlled tra- 
jectory V^^jVp^Jh? , respectively. The red dashed line at 1/4 
represents the Heisenberg uncertainty limit fi^/4. 



^{PM) satisfy 



average photon number can be estimated as: 

1 
2 

1 

" 2' 




Thus, from Eqs. (|24p and (|30p . the controlled stationary which, under the parameters given in Eqs. 



(31) 

(HZl), 
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(|28)) . and ([29)) . is smaller than the uncontrolled stationary 
average photon number 



1 



1 



Further calculations show that, if we increase Vp/ior fur- 
ther, the minimum average photon number that can be 
reached by our strategy is about 0.35 (corresponding to 
an effective temperature T^^ « 5.6 mK). 



where um is given in Eq. ([TT 

Alternatively, we can use an effective temperature T^s 
to quantify the cooling effect which is defined by: 



1 



or, equivalently, 



TcS — 



fcs In 



(32) 



The controlled stationary effective temperature T^^ can 
be estimated as follows: 



M 



off 



/cs In 



(33) 



We now give some physical interpretations of our cool- 
ing strategy. There are two competing processes that 
determine the stationary effective temperature of the 
nanomechanical beam. The cooling process is provided 
by the leakage of the transmission line resonator, whose 
energy gap is larger than ksT such that the thermal exci- 
tation from the environment could be negligible. Energy 
flows from the beam to the transmission line resonator 
via the coupling between them, and then it is dissipated 
via the leakage of the transmission line resonator. An 
opposing heating process is provided by the thermal ex- 
citation of the beam from the environment. Without 
applying quantum feedback control on the transmission 
line resonator, the cooling process of the beam is weak 
compared with the heating process, which leads to the 
failure of cooling. When we apply quantum feedback 
control and adjust the control parameters to be in the 
region given by Eq. P5|) , the decay of the beam caused 
by the leakage from the transmission line resonator is 
enhanced to overwhelm the heating process. Thus, the 
nanomechanical beam is effectively cooled. 

To show the validity of our proposal, let us show some 
numerical examples. The system parameters are chosen 
as in Eqs. ((26)) and ((27)) . and the feedback control pa- 
rameters Vx and Vp are chosen such that 



U:r = 0.5tJT, 0.3 < < 1. 



(34) 



The numerical results in Fig. [5] show that the average 
photon number and the effective temperature of the beam 
under control are reduced compared with the uncon- 
trolled case. It means that our strategy can indeed effec- 
tively cool the motion of the beam. With the parameters 
given in Eqs. P5|) . (P7)l . and ([M]) . the minimum average 
photon number that can be reached is about 0.43 (cor- 
responding to an effective temperature T^g « 6.3 mK). 




08 ^TTw^ 



FIG. 5: (color online) Cooling the nanomechanical beam 
for (a) the average photon number n, and (b) the effective 
temperature Tcff versus the normalized control parameter 
Vp = Vp/uoT- The blue solid lines representing the average 
photon number n*^ and effective temperature under con- 
trol are below the green lines with asterisks representing the 
corresponding average photon number n'"^ and effective tem- 
perature T^g without control, which means that the designed 
feedback control could effectively cool the fluctuations of the 
beam induced by thermal noises. 



VI. CONCLUSIONS 

In summary, we have investigated the possibility of us- 
ing quantum feedback control to squeeze and cool down 
the fluctuations of a nanomechanical beam embedded in a 
coupled transmission line resonator-SQUID-mechanical 
beam quantum circuit. The leakage of the electromag- 
netic field from the transmission line resonator is de- 
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tected using homodyne measurement, and the measure- 
ment output is then used to design a quantum-feedback- 
control signal to drive the electromagnetic field in the 
transmission line resonator. The designed quantum- 
feedback-control protocol indirectly affects the motion of 
the beam by the inductive coupling between the trans- 
mission line resonator and the beam via the rf- SQUID. 
After adiabatically eliminating the degrees of freedom of 
the rf-SQUID and the transmission line resonator, the 
quantum feedback control results in a two-photon term in 
the effective Hamiltonian and additional damping terms 
for the beam, which lead to squeezing and cooling for the 
beam. By varying the feedback control parameters, the 
variance of either the position or momentum of the beam 
could be squeezed, and the variance of the momentum of 
the beam could even be squeezed to be less than the stan- 
dard quantum limit h/2. Meanwhile, the average photon 
number (or, equivalently, the effective temperature) of 
the beam could be reduced effectively by applying con- 
trol, compared with the uncontrolled case. 

Although the thermal motion of the beam could 
be effectively suppressed by the proposed quantum 
feedback control protocol, our calculations show that 
the beam has not achieved the ground state. Further 
work will focus on extending our results to explore ways 
to further lower the achievable effective temperature 
or even attain the ground state of the beam. Another 
possible direction is to consider nonlinear effects [g^l 
of the nanomechanical oscillator, which may affect the 
achievable cooling temperature and the squeezing effects 
induced by the quantum feedback control. 
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APPENDIX A: DERIVATION OF THE 
QUANTUM STOCHASTIC DIFFERENTIAL 
EQUATION © 

Consider a bosonic model of the environmental degrees 
of freedom, and assume that the interactions between the 
degrees of freedom of the system and the environmental 
ones have linear interactions. Then, under the rotating- 
wave approximation, the Hamiltonian of the total sys- 
tem, composed of the degrees of freedom of the system 



and the environment, can be expressed as: 

Htot = Heft + j huJCpj^CpTduJ + J hwclrpCeTduj 

+ J hLJclj^jCeMduj + J hajclgCesduj 
+h J duj g*rp{uj)c^prpa + gpT{Lj)a^CpT 
+h [ did glj.{u})c\j,a + geT{i-o)a) CeT 
+h I duj g*j^.j{uj)clj^,jb + geM{uj)b'^CeM 



+hMr / duj 



+hM^ / duj glsi'^yls^z + ges{^)crzCeS 



where CpT (cIt), CeT (4^), CeM (c^m)' CeS {els) ^'^s, 

respectively, the annihilation (creation) operators of dif- 
ferent environmental degrees of freedom, which satisfy 

[c,{u;),c]{iu')] = S,,6{iu~oj'). 

The subscripts "pT" and "cT" represent the environmen- 
tal degrees of freedom interacting with the transmission 
line resonator being "probed" (pT) and not being probed 
(eT). The subscripts "eM" and "eS" denote the environ- 
mental (thus the "e" ) degrees of freedom interacting with 
the mechanical beam and the SQUID, respectively. 

Let A be a system operator, then the Heisenberg equa- 
tion for X can be written as: 

dX = ~l[X,Htotdt]^^[[X,Ht,tdt],Htotdt] 
n 2n^ 

= -l[X,H,s]dt + dLpTiX)+dLeT{X) 



+dLeM (X) + dL^s, dLeS, r ( A) , 



(Al) 



where 

dL^{X) 
and 

HpT — 
HeT = 



I [X, H^dt] ~^[[X, H^dt] , H^dt] 



h j dio 
h I du) 



HpM = h I dio 



glj.{uj)clj,a -f geT{i^)a^CeT 



H, 



eS,ip 



eS,r 



hM^ j duj 
hMr I duj 



9*esi^yls^=^ + 9es{(^)'^zCeS 
9*es{^)(is'^- + 9es{^^)cr+CeS 



Here, we expand dX to second-order differential terms, 
because we may meet quantum Wiener increments and 
the second-order terms cannot be omitted. 
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In order to eliminate the environmental degrees of free- 
dom corresponding to CpTito), we first solve the equation 
for Cj,t(w): 

CpT{uj) = --[cpTiuj),H tot] = -iuJCpriuj) - ig*pT{<^)a 



to obtain 

CpT(uJ,t) 



Jto 

+e-'"(*-*«)cpT(a;,io). 



(A2) 



Further, let us introduce the so-called first Markovian 
approximation to omit the frequency dependence of the 
couphng strength [ist : 



9pt{uj) 



2tt 



(A3) 



where 7pT and (ppT are independent of w. 

By substituting Eqs. and (|X3)) into dLpxiX), we 
have 



--[X,HpTdt] 



where 



du;\gpT{u;)\^ / dt'e— (*"*') 

Jto 

ia\t')[X, a]+a{t')[a\X])}dt 
-ig*pr J dwe'"(*-*°)4j,(c^, to)[X, a]dt 

+igpT[a\X] J duje-''^^*~*''^CpT{uJ,to)dt 
^(at[X, a]~a[X,a''])dt 
+ {-iV%¥e-''^-^dAl) [X,a] 
+ [a\X] (zv^e^'^^-diin) , 

dwe-*"(*-*'')cpT(t^,io) ) dt 



is the input quantum noise such that 

dA^ndll = {n(ujT) + l)dt, 

dij„diin = n{ujT)dt, 
diindiin = dA\jAl^ = 0, 



and 



n{uj) 



1 



Further, we have 



- — [[X,HpTdt\,HpTdt] = -^n{u:T)[[X,ala^]dt 



From the above analysis, it can be calculated that 
dLpTiX) = 



i 1 

- [X, HpTdt] — l^W^i Hpr] , Hpxdt] 



^{n{u;T) + l){a^[X,a] + [a^,X]a) 



IpT 



n{uT) {a[X, a)] -f [a,X]a) 



+ \/lpT[a\X]dAin^ 

where dAi„ = ie"^"'^ dA-,^. 

With the same analysis, we can calculate dLexiX), 
dLeM{X), dLes.Lp{X), and dLes,r{X). Furthermore, 
under the condition that hu;s,huJT 3> ksT, we have 
n{u!s), n{u>T) « 0. Thus, by substituting the above 
results into Eq. (|A1[) . averaging over the fluctuations 
caused by the thermal noises, and assuming that 

IT = 7pT + 7eT, 'n= ■ , 

7pT + 7eT 

we can obtain the quantum stochastic differential equa- 
tion dHl). 

In order to calculate the measurement output of the ho- 
modyne detection, let us recall that the input and output 
detection noises should be 



dAr, = 



-i= j dtJe-"^^*-'»^CpT{uJ,to)^ dt, 
-i= J duje-'"^^'-''>'>CpT{uj,h)^ dt, 



where the time to is an instant before the measurement 
commences, and the time ti is another instant after the 
measurement has finished. The measurement output is 
related to dAo^t by: 



dYt 



dAU. + e^' 



'dAo 



where (pLO is an adjustable phase introduced by the local 
oscillator of the homodync detection. From Eq. (|A2p . we 
have 

ft 



CpT{uJ,t) 



/ipTt 



-i4>pT . 



1 



V27r jt„ 
1 



dt'e 



/7pT£ 



~i4>pT _ 



dt'e 



j{t-t') 



j{t-t') 



ait') 



ait') 



V27r jt^ 
+e ^ ^>CpT{LO,ti). 
Thus, it can be calculated that 

diout - diin = ~iy/^e'"^''^a{t)dt, 
from which it can be shown that 

dYt = i^T^ate't'^^'^-'^^o) - iy/%^ae-'^'^''^-'^''°'> 
+e-'"l'"'dAl + e"'"^"dAin. 
By setting (/i^o = (ppT + 7r/2, we have 



-{n{LJT) + l)[[X, a1],a]dt. 



dYt = y^{a^ +a) + (^ie"t''-^dAin - i 
= Vvi¥{a^ +a) + {dAin + dAlj. 



f'-^dA 
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APPENDIX B: DERIVATION OF THE REDUCED 
STOCHASTIC MASTER EQUATION (USD 

Under the semiclassical approximation, we can obtain 
Maxwell-Bloch-type equations from the stochastic mas- 
ter equation (fT2| for the coupled beam-SQUID-resonator 
system (see, e.g., Ref. [6^|). Further, in the large- 
detuning regime (see Eq. ([5])), we have 



9ms 
Ams ' 



9sT 
AsT ' 



gMT ^ Wa/, UJt, Ws, 



where 



gAiT — gAisgsT 



A 



ST 



(Bl) 



(B2) 



Then, we can omit the frequency shifts of the beam, the 
rf-SQUID, and the transmission line resonator induced 
by the coupling between them. Under this condition, 
the Maxwell-Bloch-type equations for the coupled system 
obtained from the stochastic master equation can be 
expressed as: 

1 



(d,) = -ujsi^y) - 27S ( + -M: 



4 



(ay) = wsifTx) -27S (^Af^4 



d{a) 



-iujT(a)dt - gMT{o'z){b)dt ■ 



It 



'a)dt — iu{t)dt 



VlT 

h 



XtPt 



- ^\ dW, 



(B3) 



where dW has been given in Eq. ([TS 



Vx. 



{^T 



V„. 



(Pt) - {PT? 



are the variances of the normalized position and momen- 
tum operators of the transmission line resonator given by 
Eq. ini); and 



CxtPt — 



/ xtPt + Ptxt \ 
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{xt){pt) 



is the corresponding symmetric covariance. Under the 
semiclassical approximation and the condition (|B1[) . Vx^, 
Vpj. and Cxtpt can be given by the following equations: 



Vx^ = 



-JtVxt + "^^tCxtpt 



Vr 



h'jT 



Pt 



C. 



XtPt 



-ItVpt — 2u:tCxtpt "I ^ 

-2?77tC^^p^, 

-ItCxtpt + ^tV^t ~ ^tVxj. 
-2777T (Vx^ - ^ 



Gx 



(B4) 



By substituting the feedback control ((T6)) into 
Eq. (|B3[) . we can replace the last equation in (|B3p by 
the following equation: 



It 



a)dt 



d{a) = -iujT{a)dt - gMT{crz){b)dt ■ 

+ivx{a + a))dt — ivp{—ia + ia^)dt 

^ (Vx^ + iCx^p^ - ^ dW. (B5) 



If the damping rates 75 and 7t of the rf-SQUID and 
the transmission line resonator are large enough such that 



IS, It > iMTiM, 



(B6) 



we can adiabatically eliminate (58|, |59| the degrees of free- 
dom of the rf-SQUID and the transmission line resonator 
to obtain the reduced equation of the beam. In fact, in 
this case, we can obtain the following stationary variances 
from Eq. ([B4l) : 



PT 



h 

2^ ^XTPT 



from which it can be verified that the fluctuation in 
Eq. (jBSp will tend to zero. Then, in the Heisenberg pic- 
ture, one finds from the stationary solution of Eqs. (|B3p 
and (HSI that 



Cib + C2b\ 



(B7) 



where Ci, C2 are given by Eq. (^J). Substituting 
Eq. (jB7|) into Eq. we can obtain the reduced 

stochastic master equation (fT9|) for the nanomechanical 
beam. 
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